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Abstract 

The problem of large order behaviour of perturbation theory for quantum mechan- 
ical systems is considered. A new approach to it is developed. An explicit mechanism 
showing the connection between large order recursive relations and classical euclidean 
equations of motion is found. Large order asymptotics of the solution to the recursive 
relations is constructed. The developed method is applicable to the excited states, as 
well as to the ground state. Singular points of the obtained asymptotics of the per- 
turbation series for eigenf unctions and density matrices are investigated and formulas 
being valid near such points are obtained. 







1 Introduction 



A factorial growth of perturbation series coefficients which implies the divergence of 
the expansion has been studied in many papers (see, for example, 0, |3|, f|, |, |6|, 
[7], H). The corresponding large order asymptotics is usually constructed by repre- 
senting the k-th order of perturbation theory as a functional integral and applying 
the saddle-point technique to the obtained formula. Therefore, high order behaviour 
of peturbation theory is specified by saddle points expressed in terms of classical 
euclidean solutions. 

Examples of emloyment of this method are: 

i) calculating large order asymptotics for quantities like euclidean Green functions 
both in quantum mechanics and quantum field theory ||; 

ii) finding high order behaviour of energy level perturbation theory in quantum 
mechanics |5|, ||, [7|, H ; 

iii) evaluating the asymptotics of ground state eigenfunction perturbation theory 
as the number of order k tends to infinity and the wave function argument is of order 

!■ 

It is the form 

H(^gp,^gq)/g (1) 

of the dependence of the Hamiltonian on the perturbation theory parameter g that 
allows us to make use of the saddle-point method and to clarify the role of classical 
euclidean solutions. Here p is a momentum operator, q is a coordinate operator, the 
behaviour of the Hamiltonian as its argument tends to zero is like p 2 /2 + q 2 /2 (the 
harmonic oscilator Hamiltonian), so there are no terms of order 1/g, l/yfg. If the form 
of the Hamiltonian differs from ([I]), one can separate it into a leading ("classical") 
part having the form ([!]) that determines classical euclidean equations of motion and 
a correction ("quantum part") being of a form like g a f(^p, ^/gq)-, ot > 0, which can 
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influence only a pre-exponential factor (some examples are given in |7], ||). 

However, one can develop the perturbation theory not only by making use of the 
path integral representation for corresponding quantities. As a rule, when one investi- 
gates quantum mechanical case and, in particular, checks the large order asymptotics 
by numerical calculation [n|, one applies another technique, namely, one studies the 
perturbation recursive relations expressing the k-th order of perturbation theory for 
the eigenvalue and the eigenfunction in terms of preceding orders. Large orders of the 
recursive relations have not been analyzed except for the case of the quartic potential, 
TC = p 2 /2 + q 2 /2 — gq 4 , see || HQ. The explicit mechanism showing the correspondence 
between recursive relations and euclidean classical solutions has not been considered 
yet. It is studied in this paper. 

The main idea of the method is the following. First, perturbation recursive rela- 
tions is written in terms of the operator of differentiation d/dk with respect to the 
perturbation theory order k. Namely, the expression for the k-th order of the 
wave function perturbation expansion contains the preceding orders ^k-p which can 
be expressed as e _p ^\l/fc. Second, consider large order behaviour of recursive rela- 
tions, i.e. let k be proportional to some large parameter N. As it is shown in section 
2, there exists such rescaling of the wave function argument (being n-dimensional 
in the case of the n-dimensional quantum mechanics) and of the eigenfunction that 
after it the coefficients of the derivation operators become of order 1/N. This allows 
us to find the asymptotic solution to the modified equation in a form, analogous 
to the tunnel semiclassical approximation, namely, e~ NS , where the function S de- 
pending on the rescaled coordinates and on the ratio k/N satisfies the equation of a 
Hamilton- Jacobi type. As it is known, for solving such equation, one must consider 
the corresponding Hamiltonian system. As the number of arguments of the function 
S is equal to n + 1 (one more coordinate, k/N, is added), the classical Hamiltonian 
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system corresponding to recursive relations contains one more coordinate and one 
more canonically conjugated momentum. It occurs that the solutions to this (n + 1)- 
dimensional Hamiltonian system can be expressed in terms of the solutions to an 
n-dimensional Hamiltonian system associated with the euclidean equation of motion 
(see section 2 for more details). 

In order to find the solution to Hamilton- Jacobi equation, one must know bound- 
ary conditions for the function S. For finding such conditions, one can make use of 
the following observation. The k-th order of the eigenfunction perturbation theory 
can be presented as a product of a polynomial function by a Gaussian exponent (it 
is the ground state wave function of the non-perturbated system being a harmonic 
oscilator). The polynomial function can be approximated at large x by a leading 
power function. The coefficient of it is calculable explicitly. On the other hand, an 
asymptotic formula must be valid, in particular, at large values of an argument. Sec- 
tion 3 contains a comparison of these two formulas and obtaining a requirement for 
the function S, as well as for the pre-exponential factor, as their argument tends to 
infinity. This allows us to solve corresponding equation. Section 4 deals with such 
calculation. 

The obtained formulas being applicable not only to the ground state perturbation 
theory but also to the excited state perturbaton theory are valid when the eigenfunc- 
tion argument is of order y/k. When it is fixed, the asymptotics has not a saddle- 
point form ||, but, nevertheless, can be obtained by recursive relations analysis. Such 
derivation is presented in section 5. 

Section 6 contains finding the densiy matrix large order asymptotics and calcula- 
tion of some matrix elements. Section 7 deals with concluding remarks. 
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2 Perturbation recursive relations and Hamilto- 
nian systems 

The purpose of this section is to show a connection between perturbation recursive 
relations associated with the Hamiltonians having the dependence ([I]) on the pertur- 
bation theory parameter and the Hamiltonian system associated with the classical 
analog of eq. ([[]). 

Consider for the simplicity the one-dimensional quantum mechanical case, choose 
the Hamiltonian to be a sum of kinetic and potential terms and suppose the potential 
to be even, so 

n = "53? + \ V( ^ X) = "53? + T + a "*' + ^ + - (2) 

Examine the perturbation theory for the equation Tity = Ety . Both eigenvalue and 
eigenfunction can be expanded into asymptotic series, 

E n = Y,9 k E n ,k,^n{x) = J2 9** 'n,k( X ) ■ 

Here n is a number of energy level. Notice that the wave function is expanded at 
fixed value of its argument x. For the zeroth order of perturbation theory we have 
an oscilator equation, 

( "53? + T " B "' o) *"'° = a 

The n-th energy level is equal to E n ^ = n + 1/2, while the wave function has the 
form of a product of the n-th degree Hermite polynomial by the Gaussian exponent, 

* n)0 = constH n (x)e~ x2/2 . 

Let us choose the constant in this formula in order to make the coefficient of the 
leading power x n to be equal to 1, so large x behaviour of the n-th state eigenfunction 
is as follows, 

*„,o(x) ~x n e~ x2/2 ,x^ oo. (3) 



The perturbation theory can be developed by the following recursive relations ob- 
tained from eq.(||), 

(id 2 x 2 \ k 

\2d^ + YJ ^ n > k ^ + a i x4 ^n,k-l(x) + ... = J2 E ^n,k-j{x). (4) 

Notice that ^ n ,k is defined from this relation up to a factor const ^^o, so one must 
impose an additional condition on ^ n ^, for example, the following, 

d" 



dx r ' 



(e x ^ n , k )(x = 0) = 0,k>l. (5) 



Consider the construction of high order asymptotics. Suppose k to be large, i.e. to be 
proportional to some large perameter N, k = Nk. The aim is to find such rescaling of 
coordinate x and function ^ n ,k that provides a possibility of searching the asymptotics 
in a form e~ NS ^ q \ where q is a set of rescaled coordinates and k. Such rescaling of x 
is x = y/~Nr], namely, the expression (—d 2 /dx 2 + x 2 )e~ NS ^ K ' v ' ) is then approximately 
equal to N(—(dS/dr]) 2 + r] 2 )e~ NS ^ K ' v \ i.e. these two terms are of the same order. 
To make other terms in the left-hand side of eq.(^) of this order it is necessary to 
suppose that is approximately in iV times greater than *f? n) k-i, so we come to 
the following substitution, 

k = Nk, x = VNt], m n , k = N k (j) n (K, n) ~ jvV^ 5 ^, iV -> oo. (6) 

These relations show us that one can find large k asymptotics as xj \fk — > const. The 
latter condition coincides with obtained in by the path integral technique for the 
ground state. 

We can notice that all the terms of the left-hand side of recursive relations (01) play 
an important role in constructing large order asymptotics because of the dependence 
([T]) of the Hamiltonian. If there were terms such as g 2 bx A (for example, if the potential 
had a form x 2 /2 + g 2 (bx 4 + ax 6 )) then these terms would not influence a leading order 
approximation for S. This is in agreement with 0. 



One can also expect that E n ^ asymptotics has the form, 

E n>k ~ N k e- N3 <M, (7) 

where the variable r\ is substituted by zero. Namely, function \l/ n fc is presented as a 
product of a polynomial function by a Gaussian exponent, 

*n,k = J2 A n,k,lX 1 e~ X2/2 - (8) 

As the condition (|5]) is chosen, A n ^^ n = 0. It is easy to obtain from eq.(f|) that 
E n ,k = (n + l)(n + 2)A n>fcin+2 /2 + Ep= 4 a P ^n,fc,n- P , i-e. £ n ,fc is expressed as a finite 
sum (the number of terms does not depend on k) containing the derivatives of finite 
order (which also does not depend on k) of the function ^> n ^(x)e x2 l 2 as its argument 
is equal to zero. As the ty n> k asymptotics has the form @, we can make a limit rj — > 
and come to eq.(|7|). 

It will be shown that S(j], k) > 5(0, k). This implies that we can omit the terms 
in the right-hand side of eq.@ as j ~ k because they are exponentially small in 
comparison with the each order of the left-hand side. This omission is not valid only 
if q belongs to a small region near zero which decreases as k increases. 

Notice also that if j ~ 1 then the corresponding term in the right-hand side of 
eq.(^j) is of order 1/N^ +1 in comparison with all terms in the left-hand side and, 
therefore, can influence only a pre-exponential factor and corrections. Thus, the 
function S is determined only by left-hand side terms and does not depend on a 
number of energy level. 

We can also rewrite quantities such as ^ n ^~ P as e~ pd ^ 9k ^ n ,k — e~^^*^n,k an d, 
therefore, come to the following form of recursive relations, 

[H{q, -~) - ij^ol^fa, *0 = 0. ( 9 ) 

where q is a set of r\ and k, while 

1<9, Id 2 i o , r/ i__a_^ , \ 

H ^-Ny q ) = -2N^ +eN {r]e ~ 2NdK) - ( } 
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The terms of order lesser than 1/N are omitted in eq.(^). The generalization of 
formula @ to multidimensional case is straightforward. 



Let us look for the asymptotics in a form [11 



MQ) = N a X n(q)e- NS ^(l + 0(1/N)), (11) 

where a is a parameter of order 0(1). One can obtain the Hamilton- Jacobi equation 
for S, 

H(q,dS/dq) = 0. (12) 

In order to solve it, one must consider the corresponding classical Hamiltonian system. 
In the present case, the classical Hamiltonian function depends on two coordinates r\ 
and k, two momenta, p v and p K , and can be obtained from eq. ([lO|) by substituting 
the operators —N~ l d/dr] and — N~ 1 d/dK by p v and p K correspondingly, 

H(r ] ,K,p v ,p K ) = ±p 2 v -e^V( V e^ 2 ). (13) 

Let us consider it in more details. First of all, the function (|K|) does not depend on 
the coordinate k, so the momentum p K is an integral of motion. Then, the quantities 

satisfy the euclidean equations of motion, 

Q = P,P = V'(Q). (15) 

Thus, the examination of the Hamiltonian system corresponding to the perturbative 
recursive relations can be reduced to the study of the Hamiltonian system associated 
with moving of a particle in the inverted potential. If the latter Hamiltonian system 
is integrable (for example, in one- dimensional case), the Hamiltonian system fll5|) is 
integrable, too. The equation for k has the form, 



k = e~ PK 



1-9.A. 

2 dQ 



V(Q) (16) 
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and can be solved exactly. 

For finding large order asymptotics of matrix elements, it is useful to know high 
order behaviour of perturbation theory for quantities like p m,n (x,y) = ty m (x)^ n (y). 
It can be found by making use of the ^f n> k asymptotics (this technique is developed 
in section 6). But another way of constructing p m ' n (x,y) large order asymptotics is 
to apply technique developed in this section. Namely, the function p m ' n (x,y) obeys 
the equation, 

One can develop the perturbation theory p m ' n (x,y) = Y,Q k Pk ,n (x,y), change the 
variables k = Nk,x = y/Nrji,y = ^/Nr] 2 , p™' n = N k R m ^ ni rewrite the recursive 
relations corresponding to eq. (fPTD in a way analogous to eq.(H), obtain a Hamilton- 
Jacobi equation and reduce it to a Hamiltonian system that contains three coordinates 
x iViiV2, three momenta p K ,p Vl ,P V 2- The corresponding Hamiltonian function is then 

H = p 2 j2+p 2 j2 - e^V( m e^ 2 ) - e~*«V (m<? K/2 ) ■ 

The corresponding Hamiltonian equations can be solved as well as in previous case. 

One can also examine the perturbation theory for eigenfunctions in other rep- 
resentations. It occurs that perturbative recursive relations in any presentation are 
connected with the Hamiltonian system related with the investigated one by a canon- 
ical transformation. Thus, consideration of the high orders is associated with the 
investigation of the euclidean equations of motion (|i5|). 

Let us illustrate this declaration by an example. Consider for the simplicity case 
of the ground state, i.e. n = 0. The function \l/o,fc has then the form (||), 

2k 

* 0>k (x) = £ A Km x 2m e~ x2 l\ 

771=1 



s 



where the polynomial coefficients satisfy the relations: 

k 

-(m + l)(2m + l)A k , m +i + 2mA Km + a 4 A fc _i )m _ 2 + a 6 A k _ 2 , m -3 + ■■■ = Y1 E 0d A k _ j:m . 

3=1 

Assume that k and m are proportional to a large parameter, /c = Nk, m = Nfi, and 
substract the factor N k ~ m from the coefficients, A k ^ m = N k ~ m Ak,m- For the reduced 
coefficients Ak, m one can obtain in a leading order a pseudodifferential equation, 

1 9 1 d . . 

where the corresponding Hamiltonian have the form 

H(fx, K,pp,p K ) = 2(/ie- p " /2 - e p " /2 /2) 2 - e - pK V(e (p ^ )/2 ) 

and can be obtained from the investigated one by the canonical transformation: 

p„ = 2\nr],n = rj(r]-p ri )/2. 

Thus, one can apply the same technique in order to obtain high order asymptotics in 
other representations. 

3 Large x behaviour of the asymptotics and 
boundary conditions 

As it has been shown in previous section, the problem of finding the large order 
asymptotics for the wave function perturbation theory is reduced to the problem of 
solving the Hamilton- Jacobi equation flT2] ) and the equation for Xn which is obtainable 
from eq.(|Sp. To solve them, one must know boundary conditions which are to be found 
from other reasoning. 

For this purpose, consider the behaviour of ^ n ,k at large x. This allows us to 
impose boundary conditions on the functions Xn, S at large rj. 
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As ty n> k has the form @, large x behaviour of it can be obtained by approximating 
the polynomial function by the leading order, 

*njk{x) ~ c nM x ik+n e- x2 ' 2 . (18) 

Quantities c n ^ obey recursive relations, 4kc n ^ + a^c^k-i = 0,so 

c k = {-a 4 /A) k /kl 

Let us make use of eq.(|l8|) in order to find the conditions to be imposed on the 
functions S, Xn and the constant a. Changing the variables ([]) and employing the 
asymptotic formula for k\ at large k,k\ ~ (k/e) k V2Trk, we find that 

y/2TVr] 

Thus, one can expect that 

a = (n-l)/2,xn~ -f= =,S~ — -«ln( — — ) - k (19) 

V27TK 2 4k 

at larges rj. Let us discuss now when the approximation (|19"D works. In order to 
answer this question notice that the large parameter N can be chosen in different 
ways, it can be rescaled, iV — > aN. Then other variables also rescale, k —>■ n/a, 
r] — > ri/^fa, (p — > <j)a~ . As the function does not alter, we have 

Therefore, the following properties of x n ,S 

X„(k, 77) = a a Xn(K/a,r]/y/a), S(k,T]) = — rein a + aS(n/a,r]/\/a), 

which are obviously satisfied for eqs.flIDD, take place. We see that the values of these 
functions at fixed are related with each other. Thus, we notice that if the 

approximation ([19]) is valid at some value of 77, k then it is applicable at all values of 
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77, k having the same ratio T]/\/k. Therefore, tj/^/k is the only parameter determining 
whether approximation fll9| ) is valid or not. Thus, eq.([H|) is applicable at 

77/V«»l. (20) 

Taking into account eq.(^) as a boundary condition for S, Xn, find the solution to 
the Hamilton- Jacobi equation fll2f) and the equation for \n obtained also from eq.(|). 



4 Large order asymptotics for the eigenfunction 
perturbation theory 

In this section the functions S, Xn and the constant a are found, so that large order 
behaviour of the eigenfunction perturbation theory is obtained. 

4.1 Finding a solution to Hamilton- Jacobi equation 

In order to obtain a solution to the stationary Hamilton- Jacobi equation let us con- 
sider first a non- stationary equation for the function S(q, t) depending on / coordinates 
q and time parameter t, 

dS TTI dS, 

and remind a helpful property of its solution (see, for example, Consider the 

2 /-dimensional phase space with first / coordinates being momenta p and the rest 
being coordinates q. Examine the /-dimensional surface in it, 

pM = ?r ( 21 ) 

oqi 

depending on t (the arguments q,t of the function S are omitted in eq. fl2T|)). Consider 
any point (p°, q°) on the surface ([H]) when t = and its time evolution, i.e. the 
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trajectory (p , q ) in the phase space which satisfies the Hamiltonian system, 

. t _ OH , t _ OH 
Pi ~~dq~ ,(ll ~di l ' 



It can be shown |TT| that the point (p*, q l ) belongs to the surface fl21~|) at time moment 
t. 

Making use of this property let us find a solution to the stationary equation 



( |12D which obeys eq. (|19D under condition (p0|) . This means that the surface (|21|) 
transforms under this condition to the surface 



p K m( — — ),p v ~ v (l (22) 



Reconstruct now the whole surface (pij), making use of the asymptotics (|2"2"D. We 
know that the surface is invariant under time evolution, so we can consider the 
surface consisting of the trajectories that satisfy Hamiltonian systam and start from 
the asymptotical surface (^). Consider any point on it and treat it as an initial 
condition for the Hamiltonian system associated with eq. (|13D , which reduces after 
substitution flTil) to eq.(pl|). An initial condition for it is as follows, 



P ~ Q ~ 



4k 



CI477 4 



We see that if the condition (p0|) is satisfied then this point (P, Q) is closed to the 
minimum of the Hamiltonian (|l]),P = Q = 0, so that the trajectory should start from 
it. To countrbalance it, the parameter p K can take on any value under condition (PU|). 
Therefore, the surface ([H]) consists of the following treajectories, 

^t) = Q(T)e-^ 2 ,p v ( Pft ,r) = P(r)e-^ 2 , 

(23) 

K>{p K ,r) = \(T)e~ PK ,p K (p K ,T) =p K , 
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where P(t),Q(t) is a solution to eq.([15|) which starts from zero as r — > — oo, the 
parameter A(r) is found from eq.(|T6|) and condition A(— oo) = derived from eq.(^). 
It has the form 

^ )= L dT[1 'ih v(Q(T)y 

The parameter r is defined up to shifting by a constant. In order to fix it choose the 
parameter c in the r — > — oo asymptotics, Q{t) ~ ce r , to be equal to 1. 

Consider an example. Examine the function V in eq.(0) which is depicted in 
fig. la. It has a relative minimum at Q = 0, is positive as \Q\ < Q + and negative as 
\Q\ > Q + . The euclidean classical solution to eq. ([15|) starting from P = Q = as 
r — > — oo behaves as follows (fig. lb). At first stage the coordinate Q increases, then 
the solution reflects from the turning point Q + and comes back as r — > +oo to point 
P — Q — 0. The parameter A has then the form 

A(r) = s(t) - Q(t)Q(t)/2, (24) 

where s(r) is the action of the euclidean solution, s(r) = J^ OQ (Q 2 /2 + V{Q)). The 
quantity A has the meaning of the area of the section-lined region in fig. lb. The 
corresponding curves at k — r] space at various values of p K are shown in fig. 2. In 
order to convert from quantities k, rj to quantities p K , r one should first find the 
parameter r from the relation between it and £ = i]/ \/k, 



i = Q{t)/JMt) (25) 

and then obtain p K = ln(A/fi;). Notice that large values of £ correspond to r — > — oo, 
while small values of £ are associated with r — > +oo. For the simplicity, consider the 
case, when the dependence between £ and r is one-to-one, i.e. £ < everywhere. In 
general case, the solution to Hamilton- Jacobi equation is multivalued, and it is neces- 
sary to apply the Maslov canonical operator technique |Tl|] for finding the asymptotic 
solution to eq.(||). 
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Instead of the parameter r, one can use the variable Q such as Q = Q(r). The 
correspondence between £ and Q being multivalued is shown in fig.lc. 
In order to find the function S satisfying eq . consider the integral 



p K dn + p v dr] (26) 

which does not depend on the curve connecting two points and equals to 

(s(r) + \{r)p K )e"^ + const (27) 

The derivatives of it with respect to k and r) are equal to the functions p K and p v after 
expression of their arguments p K and r through the quantities K and n. As the surface 



(|23|) lies on the surface H = 0, the function (|27|) obeys the stationary Hamilton- Jacobi 
equation (O). Because of the boundary conditions JI5D, the constant in eq.(|2"7|) is 
equal to zero. Converting the parameter p K to k, r, we obtain that 

SM = K (fM + (28) 

the argument n is to be expressed through the parameter r according to eq.([25|). 
Substituting eq. (|28|) to eq. flTTD , making use of the asymptotic formula for k\ at larges 
k and applying the expression x/V~k = £, we obtain that 

*n,fc(^) ~ fc!e" fcA «), (29) 

where the pre-exponential factor is omitted, 

A(0 = ^+lnA(r)-l (30) 

Remember now that the right-hand side of the recursive relations (ffl) has been ne- 
glected. This procedure is valid if e~ kA ^ ^> e~ fcA (°)e~ fc ^ 2//2 . Consider the function 
/(0 = A(0) + £ 2 /2 - A(0- At £ = we have /(£) = 0, while /'(£) = (Q - Q)/^A. 
When Q < 0, this quantity is obviously positive. As we have supposed £ to be neg- 
ative, the quantity A should be positive at Q > 0, i.e. Q 2 — QQ > 0. This means 
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that Q < Q. Thus, the derivative /'(C) is positive everywhere, so that /(£) > 0. The 
right-hand side of eq.(|J) can be omitted when kf(£) 3> l,i.e. 

ke » i (31) 

because at small £ the function / is of order £ 2 . 

Exponential asymptotics (p9[) can be also checked by numerical calculations. Fig. 3 
contains the graphs of functions 

^) = _l ln |£^| (32 ) 

as k = 10,20,30 (dashed lines) and function (j30|) (solid line). The potential has the 
form 

V(Q)=Q 2 /2-Q\ (33) 
We see that functions A k (£) indeed tend to at larges k. 

4.2 Calculation of the pre-exponential factor 

For finding the pre-exponential factor in eq.flTTD consider its substitution to eq.(^). 
The following equation for Xn can be obtained: 



H(a 93 1 9 )+ 1 F 



X = 0, (34) 



Expand the function H into a series in 1/N ||11|| . The term of order 0(1) vanishes 
because of the Hamilton- Jacobi equation, the term of order 0(1/N) determines the 
equation for function \ n . This equation has the form: 

( 1 dH d 1 d 2 H d 2 +E \ , q s = (35) 

\ Ndp m dq m 2N dp m dp n dq m dq n n, ° J 

where the summation over repeated indices m, n is assumed, the coordinates k, r\ are 
denoted by qi,q2, the momenta p K ,p v are denoted by P\,P2- Formula (|35| ) can be 
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derived by the expansion of the function H into a series in powers of the momenta 
and by checking this formula for power functions. Let us substitute variables q for 
parameters r, p K according to (|23|) . The operator transforms then to the 

operator of differentiation with respect to r. 

In order to simplify eq. (|30| ) notice that the Jacobian of the substitution r, p K by 

D(n,K,) , , 



satisfies the equation [11 



8J_ _ d 2 H d 2 S d 2 H 
dr dp m dp n dq m dq n dp m dq m 

Making use of the nullification of the Hamiltonian second derivatives with respect 
to p m ,qm and substracting the factor J -1 / 2 from the function Xn- Xn = fn/y/~J, one 
finds the equation for the function /„ 

^ = (n + l/2)/ B . 

Therefore, the function f n , as well as Xn, is defined up to an arbitrary factor which 
depends on p K and does not depend on r. Let us use the fact that the Jacobian J 
has the following form in this case, 

3p« 



J = exp 



• Q ■ 

-AQ + |A 



2 

and substract the square root of e~ 3pK ^ 2 into a total r-independent factor. The func- 
tion Xn has then the form: 

e r(n+l/2) ( e p K ) 

Xn = , . yK =" ) (37) 



'QX/2 - XQ 

To find the function g, let us make use of the condition (18) being valid at large £, 
i.e. at t — > — oo. Under these circumstamces we have: 



Q ~ e r , V - /2 ~ -a 4 Q 4 , A ~ -a 4 e 4r /4 
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so that 



Consider now eq.(|19|) for Xn and substitute the variables r], k by r, p K according to 
eq . (|23|) . The following expression can be found, 

e (n-2)r e -^-p K 

Xr, 



Comparing two obtained expressions for Xm we find that 



n-l. 



g (e PK ) = e-— pK /V2TT. (38) 

Making use of eq.(|37|) for Xn, eq.(pHl) for the function g, the relation ( Jl9| ) for the 
constant a and eqs.(p3|), collecting all factors in the asymptotics ( |TTD and applying 
the asymptotic formula for k\, one can obtain that 

(n+l/2)r U] 

Vn&Vk) = 1 . . — 7f (fc/A)("- 1 )/ 2 e- fcA «)(l + 0(1/*))- (39) 

'QX/2- XQ 2ixyk 



An interesting feature of eq.(^) is the absence of the singularity in the turning point: 
although the velocity Q vanishes in this point, the Jacobian is finite because A ^ 0. 
Formula (^) can be also checked numerically. The graphs of the functions 

for the ground state wave function and the potential (|3~3"D , as well as the function 

M K , = .. ,4!) 

2WQA/2- XQ 



are shown in fig. 4 at k — 10, 20, 30, 40. Fig.4 shows us that the asymptotics ( |3"9"D is 
good at all values of £ except some region near the point £ = 0. The size of this region 
decreases as k increases. This is in agreement with eq.(|31"D determining the range of 
validity of the obtained asymptotics. 
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5 Large order behaviour of ^ n ^{x) at fixed x 

An interesting feature of the asymptotics (^) constructed in previous section is the 
singularity of the pre-exponential factor near the point £ = 0. Namely, small values 
of £ correspond to r — > +00. Find the behaviour of Q, A, A and the pre-exponential 
factor in eq.([39|) at these values of r. First, notice that the coordinate Q behaves 
as r — ■> +00 as follows, Q ~ ce _T . The value of the constant c can be found by 
the following procedure. Consider a small quantity go- The classical trajectory Q(r) 
(fig. lb) passes through this point, qo, at two euclidean time moments being approx- 
imately equal to lng and lnc/g because of the Q(r) asymptotics as r — > ±00. On 
the other hand, the difference between these time moments In c — 2 In q is equal to 



2/£+ dq/j2V(Q),so that 



exp2(lnQ + + [ Q+ dq{^L= - -)). (42) 



Making use of eqs.dJ,©,©, we obtain 

CP ^~ 

A ~ s(+oo) - cV 2t , £ ~ . A ~ In s(+oo) - £ 2 /2. 



's(+ooJ 

Taking into account all factors, find that at small £ the asymptotics (|3^) transforms 
to the following formula, 

e fe5 2 /2 c n+l/2 

{^Vk) n+1 2vr (s(+oo)) fc +« +1 /2' (43) 

i.e. the pre-exponential factor diverges as l/£ n+1 , where n is a number of level. 

On the other hand, each order of the eigenfunction perturbation theory is non- 
singular everywhere and, in particular, at the point £ = 0. The difficulty can be 
resolved as follows. As it has been shown in the previous section, the obtained 
asymptotic formula (|39|) may be not valid at k^ 2 ~ 1, i.e. at fixed argument of ^ H} k 
because two sides of eq.(T|) become of the same order. Therefore, the asymptotics 



under conditions 

k — > oo, x = const, (44) 
has another form than eq.(f4*3|) and is non-singular. It is investigated in this section 



and allows us to find the large order asymptotics for eigenvalues which is also studied 
in this section. 

Let us look for the asymptotics under conditions ([44]) in a form, 

c n+l/2 fc\ftn-l/2 

Take into account that at k^ 2 3> 1, fc£ 3 <C 1 both eqs.(|43l) and ( fi5|) are applicable. 
This means that the function X n behaves at large x as follows, 

X n {x) ~ e x2/2 /x n+1 , x -> oo. (46) 

Let us also look for the eigenvalue high order asymptotics in the form: 

c n+l/2 

E n>k ~ (s(+oo))fc+ „ + i/ 2 f " ( 4T ) 



and consider the substitution of eqs.(^5|),(^7|) to the recursive relations @. Notice 
that at fixed x and larges the first term at the left-hand side of eq.(|]) prevails others, 
while right-hand side contains two terms of the same order, En ^n^ and E n ^ n fl. 
Therefore, the following equation is obtained for X n : 

+ Y~ n ~ 1/2 ) Xn(x) = Sn * n >°( x ^ ^ 

It has the form of a harmonic oscilator equation with the right-hand side. 

First of all, solve eq.(f4"8|) with the boundary condition (^) at n = 0, so that 
^o,o (x) — e~ x ' 2//2 . Let us search for the solution in a form: 

X (x) = Y(x)e~ x2 / 2 , 

for Y' \x) one can obtain the boundary condition Y' (x) ~ 2e x2 and the equation 

-(Y'e- x2 )' = 2e~ x2 . 
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Thus,the solution X is the following: 

X (x) = -^- e - x2/2 f ' dx'e x ' 2 f dx"e- x "\ 

V7T JO JO 

the boundary condition ([5]) is taken into account. The quantity £ is equal to —2/^/%. 
Functions X n can be derived from the function Xq with the help of the creation 



operator x — d/dx. Namely, if the function X n satisfies eqs.(|4q),(fl6|) then the function 



X n+1 = — — {x - d/dx)X n 
n + 1 

satisfies these relations, too. The constant £ n+1 is then = 2£ n /{n + 1), because 
the normalizing factor in ^ n fl is chosen in order to satisfy the relation @). Notice 
also that we can add the factor C n ^ n f\ to the function X n . The constant C n can be 
fixed by the condition (||). Thus, we obtain 

i 2 n+1 

X n = -Ax- d/dx) n X + CA, , £ n = r-p. (49) 

n\ nl\/7i 



where 

n\ dx 



n\ dx r - 

the derivative is taken at x = 0. 

The asymptotics for eigenvalues fl4"5|) coincides with the obtained one in [^|. [7| . 

The asymptotic formula (|43| ) can be checked by numerical calculations, too. Con- 
sider the functions 

2k a 21 

X^ k {x)=Y,^e- x2 '\ (50) 

1=1 A k,l 

which are proportional to ^o t k{x) an d satisfy the condition 

X^ k {x)e x2/2 ~x 2 + ...,x^0. 

Eq. fllSj ) implies that this function tends as — > oo to the solution to eq. ( f4"6p as n = 
which is shown by dashed line in fig. 5. Functions Xq^ as k = 10, 20, 30, 40 are shown 
by solid lines for the potential ( |33"D in fig. 5. The figure shows that these functions 
indeed tend to the solution to eq. (|48|) . 
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6 Large order behaviour of perturbation theory 
for density matrix and some matrix elements 

In previous sections large order behaviour of ^ n ,k has been found. But for calculating 
high order asymptotics of matrix elements it is necessary to know such asymptotics 
for the quantities 

p nin \x 1 ,x 2 ) = ^ ni {x 1 )^ n2 (x 2 ) (51) 

Namely, matrix elements are expressed as a trace of a product of this matrix by the 
operator corresponding to the observable. When m = n 2 , the matrix (|5T|) is a density 
matrix. One can expand eq . (|5l"D into a series in g, 

oo 

p nin2 (Xl,X2) = EM in2 (^1^2) 
k=0 

and express the coefficients through the quantities ^ n ,ki 

k 

p n k in2 {x 1 ,x 2 ) = ^ num M^ ni ,k- m (x2). (52) 

This sum is estimated in subsection 6.1, while the perturbation series asymptotics for 
some matrix elements is looked for in subsection 6.2. 

6.1 Calculation of the p^ 1112 asymptotics 

Let us calculate the sum Q52] ) at large k. Since the ^ njm (x) asymptotics is already 
known as m is proportional to the large parameter N and x ~ \ N, one can estimate 
the sum ( p2[) under conditions 

k = Nk, xi = \ZNtji, x 2 = VNr] 2 - (53) 
Denoting m = Nfi and making use of eq. flTTp , we obtain: 

Nk 

PT n2 ~ E ^"-^(^r/Oxn^^-^^^e-^^ 1 ^™)) (54) 
Nn=a 
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Consider the point (1q being the point of minimum of the function 

B(fj) = S(ji,r)i) + S(k - fj,,r}2) (55) 

and suppose (jlq to belong the interval (0, k). The main contribution to the sum fl54|) 
is then given by any small interval ((Iq — b/N@,(i + b/N 13 ), where < /3 < 1/2. 
Therefore, the pre-exponential factors can be replaced by their values at point (i , 
while the function B can be substituted by the expression 

B((i) ~ BQaq) + ^B"((M Q ){(i - /i ) 2 . 

As the distance between neighbour values of (i is equal to 1/N, the sum over (i can 
be replaced by the integral iV / d(i. Making use of these observations and calculating 
the Gaussian integral over (i, we obtain that 

P*™ ~ N k+ ^^ exp(-NB y2ir/B"((i )x ni (» ,Vi)Xn 2 (K- //o,%) (56) 

If there are several minima of B with the same value, one should take all of them into 
account; expression ( |56l) should be replaced then by the sum of analogous expressions 
associated with corresponding minima. 

Let us make use of the explicit form of the functions \ n , S and simplify eq . (|56|) . 
Convert from quantities (i, 7)\ to parameters p Kl ,Tx and from quantities k — (i,r]2 
to parameters p K2 , T 2 according to eq.(|23f). The condition B (/i ) = means that 
P/ti — Pk 2 = Pk- Therefore, in order to find the p^ 1 " 2 asymptotics under conditions 
fl5"5D one must first find the parameters p K)Tl>T2 according to the following relations, 

k = (A(n) + A(r 2 ))e- pK ,^ = Q(n)e- p ^,r) 2 = Q(r 2 )e~ p ^. (57) 

The quantity B((iq) has then the form 

B = (s(n) + s(t 2 ) + (A(n) + A(r 2 )K)e- p % 
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while the second derivative of it is as follows, 

„ dp Kl (no,r)i) dp K2 (K- fi ,rj 2 ) drjx/drx dr) 2 /dT 2 
& (P-o) = ~ ~ = -jj r + 77 v ^ 



where J is the Jacobian (|36|). Collecting all obtained factors, one can find that the 



p% iri2 asymptotics under conditions (53) has the form of a product of a power function 



of N, slowly varying pre-exponential factor and rapidly varying exponent: 

p mn 2 _ N k+^1±^1 exp (_ j/VB() ) 7} (59) 

where 

e - ni+ 2 2 ~ 1 p^ e (ni+imr 1 +(n 2 +l/2)T 2 

D = g(ri)[Q(r 2 )A(r 2 )/2 - A(r a )Q(r 2 )] + Q(r 2 )[Q(r 1 )A(r 1 )/2 - A^OQfa)], 
the quantities p K ,Ti, r 2 are to be expressed from eq. 

6.2 High order asymptotics of some matrix elements 

This subsection deals with the following example of application of the obtained asymp- 



totics formula fl59|) . Large order behaviour of perturbation theory for matrix elements 
such as 

oo 

< n 2 \x mi {-d/dx) m2 \n 1 >= Y,9 k < n 2 \x mi {-dj dx) m2 \n x > k , (60) 

where n\ and n 2 are numbers of levels, is studied. It occurs that calculation of such 
asymptotics is based on the behaviour of eq.(^) at coinciding arguments, 

Vi = V2 = V- (61) 



Namely, the k-th order of perturbation theory for the quantity ( |60| ) can be rewritten 
in a form, 

' dxdy5{x - y)x mi {-d / dx) m2 p n k lU2 (x,y). (62) 
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The ?T7.2-th derivative of p with respect to x is approximately equal to the product of 
p by (^/NdBo/ dT]i) m2 . Therefore, it is necessary to know the functions -Bo, 7 only at 
the surface (|6l|). 



Let us now find the parameters n, T 2 ,p re when eq.(pl|) is satisfied. Notice first that 



one of possible choices of 77, r 2 is as follows, 17 = r 2 = r, and, therefore, 



J? /v^ = Q(r)/V2A(r). (63) 

This choice corresponds to the following selection of the parameter /x in eq.(p5[),/i = 
k/2. It is always the extremum of the function B(p). Eq.(|58|) implies that this 
extremum is minimum when Q > (fig. 6a) and maximum when Q < (fig. 6b). 
Therefore, if the parameter r obtained from eq. ( |63| ) corresponds to negative values of 
the velocity then it is necessary to find other extremum of B(p) that is minimum. 

Another possibility is as follows. Since the trajectory Q(t) passes through each 
point twice, one can consider the parameters 17, r 2 such that Q{t\) = Q(t 2 ) but 
Qiji) > 0,Q(t 2 ) < (see fig.6c),so that 



v/V^ = Q(n)/yJs(+oo). (64) 



This selection is possible when rj/y^K < Q + / J s(+oo), contrary to the previous case 



being possible as rj/y^K > Q+/ys(+oo), and corresponds to minimum of -B, because 
J{j 2 ) - J{t\) = e~ 3pK/2 Q(T 1 )s(+oo) > and, thereby, B"(/i ) > 0. Take also into 
account that we can exchange the variables /1 and k — p and obtain another minimum 
of B. 

Making use of the cited observations, let us find the constants B,j in the 
p™ 2 {\/Nri,v r Nrj) asymptotics: 



(a) rj/yfc > Q + /Js(+oo) : 



'A(r) k 
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(ni+ri2-l)/2 e (ni+n2+l)r 

7 



K y 
V2A7 



^2Q(r)[Q(r)A(r)/2-A(r)Q(r)]' 
Notice that in this case 5 = -QJ/X 2 < 0, so that as r decreases, -B increases. 
Therefore, the main contribution to the integral ( |62|) is given by a region near the 



critical value of rj/^/n, Q + / Ws(+oo) 



(b) 77/ < Q+/yJs(+oo) : 

5 = «(l + ln^ 



^ e (ni+l/2)ri+(n 2 +l/2)r 2 _|_ e (ni+l/2)r 2 +(n2+l/2)ri j j-gc^ 

27rs(+oo)g(ri) 

the existence of two minima of the function B is taken into account. One can notice 
that the region (a) gives an exponentially small contribution to eq.(|62|). 

Find now large order asymptotics of the matrix element (^). Its evaluation is 
based on the following observations: 

(i) one can change the parameter x by a parameter r 1; according to eq.(|48|), 
the measure dx transforms then to \/NQ(ri)e~ PK ' 2 dri, while the derivation operator 
—d/dx should be substituted by the operator of multiplication by viV P(Ti)e~ PK ^ 2 at 
the first term in eq.(|65|) and by \fNP{T2)e~ PK ^ 2 at the second; 

(ii) the following relations are satisfied: Ne~ PK = fc/s(+oo), T\ + r 2 = lnc, where 
the parameter c has the form (fl2j); 

(iii) one can divide the obtained integral over n, where T\ < t ,Q(t ) = Q + into 
two parts according to eq.(^), denote by r the quantity T\ in the first integral and r 2 
in the second one, the rest parameter is then equal to lnc — r; the integrands become 
coinciding then, while the region of integration is r < tq in the first integral and 
t > To in the second, so their sum is an integral over —00 < r < +00; 

(iv) one can employ the asymptotic formula for k\; 

(v) notice that at odd m x + m,2 + n\ + n 2 the quantity ( |60D is obviously equal to 
zero, while at even m\ + m 2 + n\ + n 2 the integrals over regions x < and x > 
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coincides, so one should increase the value of the integral over x > in 2 times and 
obtain the following asymptotic formula: 

n l + n 2 ~l~ m l "l" m 2 

< n 2 \x m \-d/dx)^\ ni > fe ~ ^ } 

7r(s(+oo)) fc ys(+oo) y 

/+00 
rfre (n2+1/2)lnc Q mi (r)P m2 (r)e (ni ~ n2)r , (66) 
-oo 

Let us discuss the obtained result (|66"1). Notice first that if we substitute r by lnc — 
r, this expression will be multiplied by (— l) m2 , while the numbers n\ and n 2 will 
exchange. This fact is in agreement with the definition flBOp. 

Second, notice that when m\ and m 2 increase, the growth of the coefficients is 
faster. The (k — l)-th order for the matrix element of x m and the fc-th order for such 
quantity for x m ~ 2 have the same order of growth. Therefore, the k-th order behaviour 
for all powers of gx 2 , g(d/dx) 2 is of the same order. 

Large order asymptotics for quantities such as " Green functions" 

< n 2 \x(r 1 , ...,r m \ni > k , (67) 

where x{t) is an operator of the form x{t) = e HT xe~ HT can be also studied. Namely, 
denote the solution to eq. (|H|) satisfying initial conditions -P(O) = P ,Q(0) = Q by 
P{ T , Qoy Pq): Q( r ) Qo, Pq)- The operator x{r) is then approximately equal to 

x{t) ~ g~ 1/2 Q(r, Jgx, -y/gd/dx). 

One can expand this function into a series in operators x and —d/dx and apply 
eq. (|66|) . It occurs then that large order asymptotics for eq.(|67D can be obtained 
from eq.(|66D with the help of substituting m x by m,m 2 by 0, Q m {j) by a product 
Q{t + ti)...Q(t + r m ).This result is in agreement with the Lipatov method of 
finding asymptotics for such quantities as (^) for the ground state, namely, eq.( 
can be presented as an integral over trajectories 
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after substitution x^fg = Q, g = X/k this integral becomes of the saddle-point type, 
saddle points being Q = Q{r + a), X = s(+oo) are to be taken into account, so the 
asymptotics of eq.(|67|) is expressed through the integral / daQ{r\ + a)...Q(r n + a). 



7 Conclusions 

In this paper the explicit mechanism showing the connection between large order re- 
cursive relations for stationary perturbation theory and classical euclidean equations 
of motion is found and high order asymptotics of the solution to the recursive relations 
is constructed. The results obtained by the path integral approach are indeed repro- 
duced in the developed method. Moreover, when one finds high order asymptotics 
of matrix element (pOD, one can also find what region of integration in eq.(^) and 
what terms in the sum ( [52]) give the main contribution to the quantity (|60|). Making 
use of the obtained large order asymptotics for the wave function and density matrix, 
one can also find the main values of quantities like e- k f( x /vk) a ^ ^he fc-th order of 
perturbation theory as k — > oo. All calculations are valid both for the ground and 
for the excited states. 

An interesting feature of the eigenfunction large order perturbation theory as 
x/Vk = const is the divergence of the pre-exponential factor at small £. This difficulty 
is analogous with the divergence of the pre-exponential factor in the semiclassical 
expansion large order asymptotics discussed in [O]. The resolve of the difficulty is 



the following. At small £ associated with the values of parameter x of order 0(1) the 



asymptotics is written in another form ([45]). It is expressed through the growing at 



infinity solution to the harmonic oscilator equation with non-trivial right-hand side 
fl48D . The connection between these asymptotics allows us to find the eigenvalue large 
order behaviour coinciding with previous papers |7j]. 

The author is indebted to V.A.Rubakov for helpful discussions. The work is 
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Figure captions 
Fig.l. 

(a) : The potential V versus Q; 

(b) : classical euclidean solution in phase space; the area of the section-lined region 
is equal to the quantity A, eq. (|24l) ; 

(c) : the correspondence between the parameter £, eq.(^), and the coordinate Q. 
Fig.2. 

A set of classical trajectories in k — r] space. 
Fig.3. 

The functions v4fc(0 defined by eq. (|32|) at various numbers of the order of pertur- 
bation theoryfc = 10, 20, 30, are shown by dashed lines; the function A(£), eq. (|30|) , is 
shown by solid line. 

Fig.4. 



The ratio M k (£), eq.©, at k = 10, 20, 30, 40 and the function M(£) specified by 

eq-©- 
Fig.5. 

The functions X 0j k(x) defined by eq.(|50D at k — 10, 20, 30, 40 (solid lines) and the 
solution to eq.([48D (dashed line). 
Fig.6. 

Classical trajectories in k — rj space which correspond to various choices of the 
parameter //, an extremum point of the function fl55|) : 

(a) fi = k/2, Q(r) > 0; 

(b) n = k/2, Q(t) < 0; 

(c) /X ? k/2. 
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